Incorporating chemical signalling factors into cell-based models of growing epithelial tissues

In this paper we present a comprehensive computational framework within which the effects of chemical signalling factors on growing epithelial tissues can be studied. The method incorporates a vertex-based cell model, in conjunction with a solver for the governing chemical equations. The vertex model provides a natural mesh for the finite element method (FEM), with node movements determined by force laws. The arbitrary Lagrangian–Eulerian formulation is adopted to account for domain movement between iterations. The effects of cell proliferation and junctional rearrangements on the mesh are also examined. By implementing refinements of the mesh we show that the finite element (FE) approximation converges towards an accurate numerical solution. The potential utility of the system is demonstrated in the context of Decapentaplegic (Dpp), a morphogen which plays a crucial role in development of the Drosophila imaginal wing disc. Despite the presence of a Dpp gradient, growth is uniform across the wing disc. We make the growth rate of cells dependent on Dpp concentration and show that the number of proliferation events increases in regions of high concentration. This allows hypotheses regarding mechanisms of growth control to be rigorously tested. The method we describe may be adapted to a range of potential application areas, and to other cell-based models with designated node movements, to accurately probe the role of morphogens in epithelial tissues

Biophysically Realistic Filament Bending Dynamics in Agent-Based Biological Simulation

An appealing tool for study of the complex biological behaviors that can emerge from networks of simple molecular interactions is an agent-based, computational simulation that explicitly tracks small-scale local interactions – following thousands to millions of states through time. For many critical cell processes (e.g. cytokinetic furrow specification, nuclear centration, cytokinesis), the flexible nature of cytoskeletal filaments is likely to be critical. Any computer model that hopes to explain the complex emergent behaviors in these processes therefore needs to encode filament flexibility in a realistic manner. Here I present a numerically convenient and biophysically realistic method for modeling cytoskeletal filament flexibility in silico. Each cytoskeletal filament is represented by a series of rigid segments linked end-to-end in series with a variable attachment point for the translational elastic element. This connection scheme allows an empirically tuning, for a wide range of segment sizes, viscosities, and time-steps, that endows any filament species with the experimentally observed (or theoretically expected) static force deflection, relaxation time-constant, and thermal writhing motions. I additionally employ a unique pair of elastic elements – one representing the axial and the other the bending rigidity– that formulate the restoring force in terms of single time-step constraint resolution. This method is highly local –adjacent rigid segments of a filament only interact with one another through constraint forces—and is thus well-suited to simulations in which arbitrary additional forces (e.g. those representing interactions of a filament with other bodies or cross-links / entanglements between filaments) may be present. Implementation in code is straightforward; Java source code is available at www.celldynamics.org

Microdevices for extensional rheometry of low viscosity elastic liquids : a review

Extensional flows and the underlying stability/instability mechanisms are of extreme relevance to the efficient operation of inkjet printing, coating processes and drug delivery systems, as well as for the generation of micro droplets. The development of an extensional rheometer to characterize the extensional properties of low viscosity fluids has therefore stimulated great interest of researchers, particularly in the last decade. Microfluidics has proven to be an extraordinary working platform and different configurations of potential extensional microrheometers have been proposed. In this review, we present an overview of several successful designs, together with a critical assessment of their capabilities and limitations

A Secure Semi-Field System for the Study of Aedes aegypti

Novel vector control strategies require validation in the field before they can be widely accepted. Semi-field system (SFS) containment facilities are an intermediate step between laboratory and field trials that offer a safe, controlled environment that replicates field conditions. We developed a SFS laboratory and cage complex that simulates an urban house and yard, which is the primary habitat for Aedes aegypti, the mosquito vector of dengue in Cairns Australia. The SFS consists of a Quarantine Insectary Level-2 (QIC-2) laboratory, containing 3 constant temperature rooms, that is connected to two QIS-2 cages for housing released mosquitoes. Each cage contains the understory of a “Queenslander” timber house and associated yard. An automated air conditioning system keeps temperature and humidity to within 1°C and 5% RH of ambient conditions, respectively. Survival of released A. aegypti was high, especially for females. We are currently using the SFS to investigate the invasion of strains of Wolbachia within populations of A. aegypti

Intracellular Spatial Localization Regulated by the Microtubule Network

The commonly recognized mechanisms for spatial regulation inside the cell are membrane-bounded compartmentalization and biochemical association with subcellular organelles. We use computational modeling to investigate another spatial regulation mechanism mediated by the microtubule network in the cell. Our results demonstrate that the mitotic spindle can impose strong sequestration and concentration effects on molecules with binding affinity for microtubules, especially dynein-directed cargoes. The model can recapitulate the essence of three experimental observations on distinct microtubule network morphologies: the sequestration of germ plasm components by the mitotic spindles in the Drosophila syncytial embryo, the asymmetric cell division initiated by the time delay in centrosome maturation in the Drosophila neuroblast, and the diffusional block between neighboring energids in the Drosophila syncytial embryo. Our model thus suggests that the cell cycle-dependent changes in the microtubule network are critical for achieving different spatial regulation effects. The microtubule network provides a spatially extensive docking platform for molecules and gives rise to a “structured cytoplasm”, in contrast to a free and fluid environment

Performance of CMS muon reconstruction in pp collision events at sqrt(s) = 7 TeV

The performance of muon reconstruction, identification, and triggering in CMS has been studied using 40 inverse picobarns of data collected in pp collisions at sqrt(s) = 7 TeV at the LHC in 2010. A few benchmark sets of selection criteria covering a wide range of physics analysis needs have been examined. For all considered selections, the efficiency to reconstruct and identify a muon with a transverse momentum pT larger than a few GeV is above 95% over the whole region of pseudorapidity covered by the CMS muon system, abs(eta) < 2.4, while the probability to misidentify a hadron as a muon is well below 1%. The efficiency to trigger on single muons with pT above a few GeV is higher than 90% over the full eta range, and typically substantially better. The overall momentum scale is measured to a precision of 0.2% with muons from Z decays. The transverse momentum resolution varies from 1% to 6% depending on pseudorapidity for muons with pT below 100 GeV and, using cosmic rays, it is shown to be better than 10% in the central region up to pT = 1 TeV. Observed distributions of all quantities are well reproduced by the Monte Carlo simulation.Comment: Replaced with published version. Added journal reference and DO

Performance of CMS muon reconstruction in pp collision events at sqrt(s) = 7 TeV

The performance of muon reconstruction, identification, and triggering in CMS has been studied using 40 inverse picobarns of data collected in pp collisions at sqrt(s) = 7 TeV at the LHC in 2010. A few benchmark sets of selection criteria covering a wide range of physics analysis needs have been examined. For all considered selections, the efficiency to reconstruct and identify a muon with a transverse momentum pT larger than a few GeV is above 95% over the whole region of pseudorapidity covered by the CMS muon system, abs(eta) < 2.4, while the probability to misidentify a hadron as a muon is well below 1%. The efficiency to trigger on single muons with pT above a few GeV is higher than 90% over the full eta range, and typically substantially better. The overall momentum scale is measured to a precision of 0.2% with muons from Z decays. The transverse momentum resolution varies from 1% to 6% depending on pseudorapidity for muons with pT below 100 GeV and, using cosmic rays, it is shown to be better than 10% in the central region up to pT = 1 TeV. Observed distributions of all quantities are well reproduced by the Monte Carlo simulation.Comment: Replaced with published version. Added journal reference and DO

Azimuthal anisotropy of charged particles at high transverse momenta in PbPb collisions at sqrt(s[NN]) = 2.76 TeV

The azimuthal anisotropy of charged particles in PbPb collisions at nucleon-nucleon center-of-mass energy of 2.76 TeV is measured with the CMS detector at the LHC over an extended transverse momentum (pt) range up to approximately 60 GeV. The data cover both the low-pt region associated with hydrodynamic flow phenomena and the high-pt region where the anisotropies may reflect the path-length dependence of parton energy loss in the created medium. The anisotropy parameter (v2) of the particles is extracted by correlating charged tracks with respect to the event-plane reconstructed by using the energy deposited in forward-angle calorimeters. For the six bins of collision centrality studied, spanning the range of 0-60% most-central events, the observed v2 values are found to first increase with pt, reaching a maximum around pt = 3 GeV, and then to gradually decrease to almost zero, with the decline persisting up to at least pt = 40 GeV over the full centrality range measured.Comment: Replaced with published version. Added journal reference and DO

X-ray emission from the Sombrero galaxy: discrete sources

We present a study of discrete X-ray sources in and around the bulge-dominated, massive Sa galaxy, Sombrero (M104), based on new and archival Chandra observations with a total exposure of ~200 ks. With a detection limit of L_X = 1E37 erg/s and a field of view covering a galactocentric radius of ~30 kpc (11.5 arcminute), 383 sources are detected. Cross-correlation with Spitler et al.'s catalogue of Sombrero globular clusters (GCs) identified from HST/ACS observations reveals 41 X-rays sources in GCs, presumably low-mass X-ray binaries (LMXBs). We quantify the differential luminosity functions (LFs) for both the detected GC and field LMXBs, whose power-low indices (~1.1 for the GC-LF and ~1.6 for field-LF) are consistent with previous studies for elliptical galaxies. With precise sky positions of the GCs without a detected X-ray source, we further quantify, through a fluctuation analysis, the GC LF at fainter luminosities down to 1E35 erg/s. The derived index rules out a faint-end slope flatter than 1.1 at a 2 sigma significance, contrary to recent findings in several elliptical galaxies and the bulge of M31. On the other hand, the 2-6 keV unresolved emission places a tight constraint on the field LF, implying a flattened index of ~1.0 below 1E37 erg/s. We also detect 101 sources in the halo of Sombrero. The presence of these sources cannot be interpreted as galactic LMXBs whose spatial distribution empirically follows the starlight. Their number is also higher than the expected number of cosmic AGNs (52+/-11 [1 sigma]) whose surface density is constrained by deep X-ray surveys. We suggest that either the cosmic X-ray background is unusually high in the direction of Sombrero, or a distinct population of X-ray sources is present in the halo of Sombrero.Comment: 11 figures, 5 tables, ApJ in pres

Compressed representation of a partially defined integer function over multiple arguments

In OLAP (OnLine Analitical Processing) data are analysed in an n-dimensional cube. The cube may be represented as a partially defined function over n arguments. Considering that often the function is not defined everywhere, we ask: is there a known way of representing the function or the points in which it is defined, in a more compact manner than the trivial one